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erating universe 
stars may exist 



I. INTRODUCTION 

As alternatives to black holes, gravastars have received some attention recently par- 
tially due to the tight connection between the cosmological constant and a currently accel- 
, although very strict observational constraints on the existence of such 
In the original model of Mazur and Mottola (MM) [4;], gravastars consist 
of five layers: an internal core < r < ri, described by the de Sitter universe, an interme- 
diate thin layer of stiff fluid ri < r < r2, and an external region r > r2, described by the 
Schwarzschild solution, 

ds^ = -f(r)dt^ + T^rfr^ + (dO^ + sin^ Odif^) , (1.1) 
f{r) 

where the function /(r) is given by /(r) = 1 — 2Ai/r, in the units where c = 1 = G. 
In addition, in such a setup, two infinitely thin shells also appear, respectively, on the 
hypersurfaces r = ri and r = r2. By properly choosing the free parameters involved, one 
can show that the two shells can have only tensions but with opposite signs jj]. Visser 
and Wiltshire (VW) argued that such five-layer models can be simplified to three-layer ones 
[5I, in which the two infinitely thin shells and the intermediate region are replaced by one 
infinitely thin shell, so that the function /(r) in the metric (11.11) is given by 

1_2M r>a(r), 
m ={ (1.2) 
1 - (f) , r< a(r). 



where r = a(r) is a timelike hypersurface, at which the infinitely thin shell is located, and 
r denotes the proper time of the thin shell. The constant / = a/3/A denotes the de Sitter 
radius. On the hypersurface r = aij) Israel junction conditions yield 

id2 + \/(a) = 0, (1.3) 

where an overdot denotes the derivative with respect to the proper time r of the thin shell. 
Therefore, in the region r > a(r) the spacetime is locally Schwarzschild, while in the region 
r < a(r) it is locally de Sitter. These two different regions are connected through a dynamical 
infinitely thin shell located at r = a(r) to form a new spacetime of gravastar. 

Two different types of stable gravastars are identified by VW, stable gravastars and 
"bounded excursion" gravastars. 
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Stable gravastars: In this case, there must exist a radius Oq such that 

V{ao) = 0, V'{ao)=0, V" (ao) > 0, (1.4) 

where a prime denotes the ordinary differentiation with respect to the indicated argument. 
If and only if there exists such an for which the above conditions are satisfied, the model 
is said to be stable. Among other things, VW found that there are many equations of state 
for which the gravastar configurations are stable, while others are not 5|. Carter studied 

n 

the same problem and found new equations of state for which the gravastar is stable [6[, 
while De Benedictis et al j3] and Chirenti and RezzoUa investigated the stability of the 
original model of Mazur and Mottola against axial-perturbations, and found that gravastars 
are stable to these perturbations. Chirenti and RezzoUa also showed that their quasi-normal 
modes differ from those of a black hole of the same mass, and thus can be used to discern 
a gravastar from a black hole. Early work on dynamical thin shells with de Sitter interior 
and (A)dS-RN exterior can be found in [91]. 

"Bounded excursion" gravastars: As VW noticed, there is a less stringent notion of 
stabihty, the so-called "bounded excursion" models, in which there exist two radii ai and 02 
such that 

y(ai) = 0, r(ai)<0, V^(a2)=0, (02) > 0, (1.5) 

with V{a) < for a G (oi, 02), where 02 > a^. 

Lately, we studied this type of gravastars [l^, and found that, among other things, such 
configurations can indeed be constructed, although the region for the formation of this type 
of gravastars is very small in comparison to that of black holes. 

In this paper, we generalize our previous work [lo| to the case where the equation of state 
of the infinitely thin shell is given by p = (1 — 7)0" with 7 being a constant. When 7 = it 



reduces to the case we studied in [lO| . We shall first construct three- layer dynamical models, 
in analogy to the VW models, and then show both stable gravastars of the both types and 
black holes exist for 7 < 1 and A 7^ 0. However, when 7 > 1 even with A 7^ 0, only black 
holes are found. In the phase space, the region of gravastars and the one of black holes 
are non-zero, although the former is much smaller than the latter. The rest of the paper is 
organized as follows: In Sec. II we shall study various cases, in which all the possibilities 
of forming black holes, gravastars, de Sitter, and Minkowski spacetime exist. In Sec. Ill we 
present our main conclusions. 
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II. FORMATION OF GRAVASTARS AND BLACK HOLES FROM 
GRAVITATIONAL COLLAPSE OF PROTOTYPE GRAVASTARS 

For spacetimes given by, 

c+ {—c+f+dv+ + 2dr+) dv+ + r^d^fi, r > a(r), 



Lake found that the Israel's junction conditions yield two independent equations [ll| 



i'-i^nu-f-f-liu^f-)^i^j-, (2.2) 

M + 87radp = 47ra2[T^^M"n^], (2.3) 

where M = Aira'^a, a denotes the energy density of the shell, and p its pressure. In this 
work, we shall consider the case where inside the shell the spacetime is de Sitter, and outside 
it is Schwarzschild, namely 

/+W = 1-^, /-W = l-(y) , (2.4) 

as that studied in [lol]. The only difference is the equation of state of the thin shell, where 
in the present paper we consider a more general case in which it is taken as, 

P={l-l)a, (2.5) 



with 7 being a constant. When 7 = it reduces to the special case studied in [10|. Since 
we find 

[T.^u^n''] = TXu+^n^" - T^^u-^n-'' = 0. (2.7) 

Then, from (12. 3p we find 

- = -2(2-7)-, (2.8) 
a a 



which has the general solution, 

•'^On 2(2-7) 

^ — >^uv 

Setting 



a = aoi^r^'-'l (2.9) 
a 



M = mLo, a(r) = R{t)Lo, I = LLq, (2.10) 
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where 

Lo^(4vraoaf-))"^, (2.11) 
we find that Eg. (12.21) can be cast in the form, 

^R*^ + V{R,m,L,-f) = 0, (2.12) 

where R* = dR/d{L^^T), and 

V(R,rn,L,-<) = |-1 + - + + + - + ^| . 

(2,13) 

Therefore, for any given constants m, L and 7, Eq. (l2.12i ) uniquely determines the collapse 
of the prototype gravastar. Depending on the initial value Rq, the collapse can form either a 
black hole, or gravastar, or a Minkowski, or a de Sitter space. In the last case, the thin shell 
first collapses to a finite non-zero minimal radius and then expands to infinity. To guarantee 
that initially the spacetime does not have any kind of horizons, cosmological or event, we 
must restrict Rq to the range, 

2m<Ro< L, (2.14) 

correspondingly ao G {2A4,l). When m = = A, the thin shell disappears, and the whole 
spacetime is Minkowski. So, in the following we shall not consider this case. 
From Eq. (l2.13p . we find that 

OR 2 J ^ 2L2 

+ ^(3-27)i?^^-^. (2.15) 

As both equations, = and dV{R, m, L, 'y)/dR = 0, are quadratic in m, we can easily 
find m from these two equations, which is given by 

m,(i?,L,7) = —— — ^- -(-8R^^+^L^ + 8R^^+^L^-f 

^ ' ' 2L2(-5L2i?47 + 4L27i?47_3i?8) V ' 

-3 i?" - 4 R^^-^L^-f + 5 R^^-^L^ 

^2i?47+3^2 _4^47+3^2^^^ (\/ = V' = 0). (2.16) 

Substituting this expression in y = we find six functions L{R,'y). Due to the complexity 
of these expressions we shall not give them here explicitly. Instead, in the following we 
consider some particular cases. 
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A. m = 



In this case, the spacetime outside the thin shell is flat, and the mass of the shell com- 
pletely screens the mass of the internal de Sitter spacetime. From Eg. (12. 131) we flnd that 

Then, V'{R) = yields, 

(27 - 3) + 2i?^(2-7)^2 _ _ 5) /?8(2-7) = 0, (2.18) 

which has real solution only for 7 < 3/2 or 7 > 5/2, and the corresponding solution is given 

by 

7^[3/2,5/2). (2.19) 

Substituting the above expression into the equation V{R) = 0, we find 

1 

3-27 

(2.20) 

Figs. 1 and 2, and TableUshow the functions of Rd'j) and Ld^j), where for 7 G [3/2, 5/2), 
the equations V{R, L, 7) = and V'{R, L,'~f) =0 have no real solutions. In particular, when 
7 = 0.5 we find that ^ 0.9185586537 and R^ ^ 0.8660254039. For L < the potential 
is strictly negative as shown in Fig. 3. Thus, if the star starts to collapse at i? = Rq, it will 
collapse continuously until R = 0, whereby a Minkowski spacetime is formed. When L = Lc, 
since Rq < Lc, we can see that, the star will collapse until the center and turns the whole 
spacetime into a Minkowski. For L > Lc, the potential V{R) is positive between Ri and 
i?2, where Ri^2 are the two real roots of the equation V{R, L > Lc) = with R2 > Ri > 0. 
In this case, if the star starts to collapse with Rq < Ri, as can be seen from Fig. 3, it will 
collapse to -R = whereby a Minkowski spacetime is finally formed. If it starts to collapse 
with Ro > R2, it will first collapse to R = R2 and then starts to expand until R = 00, and 
the whole spacetime is finally de Sitter. 

It should be noted that in this case we always have V"{Rc, Lc,'j) < 0, that is, no stable 
stars exist. 



27-5 



27-3 



Rr 



2(7 - 2) 



27-5 
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TABLE I: Some values of Rc and Lc as a function of 7 for A 7^ and m = 0. 



7 
/ 


Rc 


Lc 


0.0 


0.9283177667 


0.9587624669 


0.1 


0.9199519160 


0.9535623390 


0.2 


0.9100298584 


0.9473109395 


0.3 


0.8981404718 


0.9397043233 


0.4 


0.8837277541 


0.9303202778 


0.5 


0.8660254039 


0.9185586537 


0.6 


0.8439546889 


0.9035435931 


0.7 


0.8159607706 


0.8839575014 


0.8 


0.7797433156 


0.8577474304 


0.9 


0.7318025503 


0.8215821352 


1.0 


0.6666666667 


0.7698003592 


1.1 


0.5756295655 


0.6922956261 


1.2 


0.4452231361 


0.5703420705 


1.3 


0.2598914458 


0.3713508347 


1.4 


0.0482828420 0.08734694771 


2.6 


0.4428889922 


0.8012163331 


2.7 


0.5933418502 


0.8478077860 


2.8 


0.6857500883 


0.8784631649 


2.9 


0.7485495080 


0.9002622198 


3.0 


0.7937005260 


0.9164864251 


5.0 


0.9742903329 


0.9881107444 




B. A 


= 



In this case, Eq.f l2.13p reduces to 

V{R, m, 7) = (4 i?^ - 4 - i?^^ - 4 , (2.21) 
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FIG. 1: The function Rc defined in Eq.l^JEi) for m = 0. 
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FIG. 2: The function Lc defined in Eq. ()2.19p for m = 0, where Lr^ 

from which we find that the equations V{R) = and V'{R) = have the exphcit solutions, 

1 

4-47 



5-47 



(2 7-3) 



rrir 



2(47-5) 



+ 



4 



-K 



(2.22) 
(2.23) 



5-47 

Figs. 4 and 5 show the dependence of Rc and rric on 7. The results of our previous work 
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FIG. 3: The potential V(R) for m = and 7 = 0.5. The top hne is for L > Lc, the middle line is 
for L = Lc, and the bottom line is for L < Lc- 
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FIG. 4: The function Rc defined in Eq.^Mi) for A = 0. 

are obtained when 7 = [lo|. Some representative cases are 7 = 0.5 and 7 = 3.0. For 
7 = 0.5, we find that ^ 0.5443310540 and i?c ~ 1.224744871, and for 7 = 3.0 we find 
that nic ~ 0.5120894280 and Rc ~ 1.045515917. In both cases, for m > tUc the potential 
V{R) is strictly negative as shown in Figs. 6 and 7. Then, the collapse always forms black 
holes. For m = rric, there are two different possibilities, depending on the choice of the 
initial radius Ro- In particular, if the star begins to collapse with Rq > Rc, the collapse will 
asymptotically approach the minimal radius Rc- Once it collapses to this point, the shell 
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1/2 



1 5/4 3/2 



FIG. 5: The function defined in Eq.jS^S]) for A = 0. 




FIG. 6: The potential V(R) for 7 = 0.5. The top hne is for m < nic, the middle line is for m = rric, 
and the bottom line is for m > rric. 



will stop collapsing and remains there for ever. However, in this case this point is unstable 
and any small perturbations will lead the star either to expand for ever and leave behind 
a flat spacetime, or to collapse until i? = 0, whereby a Schwarzschild black hole is finally 
formed. On the other hand, if the star begins to collapse with 2mc < Rq < Rc as shown in 
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FIG. 7: The potential V(R) for 7 = 3.0. The top Hne is for m < rric, the middle line is for m = rric, 
and the bottom line is for m > rric- 



Figs. 6 and 7, the star will collapse until a black hole is formed. For m < ttIc, the potentials 
V{R) for each case have a positive maximal, and the equation V{R,'m < rric) — has two 
positive roots Ri^2 with R2 > Ri > 0. There are two possibihties here, depending on the 
choice of the initial radius Rq- If Ro > R2, the star will first collapse to its minimal radius 
R = R2 and then expand to infinity, whereby a Minkowski spacetime is finally formed. If 
2m < Rq < Ri, the star will collapse continuously until R = 0, and a black hole will be 
finally formed. 

It should be noted that, similar to the last case, now we always have V"{Rc, rric, l) < 0, 
which means that no stable stars exist in this case, too. 



C. m / and A / 

As mentioned before, the analytic expression for Lc in the present case are too complicated 
to write out here. Instead, we shall study it numerically. Our main strategy is to start with 
rUc obtained for the case A = 0, and then gradually turn on A. We plot the potential 
V{R,mc{'y), L,j) as a function of R for any given 7, by finely tuning L until a stable 
gravastar or a bounded excursion gravastar is found [see Figs. 8-11]. The value as shown 
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TABLE II: Some values of Rc and rric as a function of 7 for A = and m 7^ 



I 


Rc 


rric 


0.0 


1.077217345 


0.5170643255 


0.1 


1.091490000 


0.5199536485 


0.2 


1.110255191 


0.5236577775 


0.3 


1.135692224 


0.5285215890 


0.4 


1.171542463 


0.5350989795 


0.5 


1.224744871 


0.5443310540 


0.6 


1.309606309 


0.5579387820 


0.7 


1.460581840 


0.5794043665 


0.8 


1.784674184 


0.6169244095 


0.9 


2.840468889 


0.6956250340 


1.0 


00 


1.000000000 


1.1 


1.475575893 


1.660022879 


1.2 


0.09921256570 


1.190550789 


1.3 


0.01134023029 


0.1360827636 


1.4 


0.007415771480 


0.006591796870 


1.5 


Undefined 


Undefined 


1.6 


3.625777848 


14.80525970 


1.7 


3.017962339 


1.043246245 


1.8 


1.867291683 


0.7407438105 


1.9 


1.501970633 


0.6398928165 


2.0 


1.333333333 


0.5925925920 


3.0 


1.045515917 


0.5120894280 


4.0 


1.017554577 


0.5045725165 


1.0x10^ 


1.000000000 


0.5000000000 



in Table III, is obtained numerically when a stable gravastar is found for a specific pair 
(^c,7)- 

It can be shown that both types of gravastars can be formed for 76 [0, 1). But for 7 > 1, 
we find that for any given values of L and m only black holes can be formed , for example, 
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TABLE III: Some values of rric and Lc obtained numerically as a function of 7 in the general case 
where A 7^ and m / 



1 


nic 




0.0 


0.5170643255 


2.8743398 


0.1 


0.5199536485 


3.1000618 


0.2 


0.5236577775 


3.3917341 


0.3 


0.5285215890 


3.7828681 


0.4 


0.5350989795 


4.3336020 


0.5 


0.5443310540 


5.1626297 


0.6 


0.5579387820 


6.5372013 


0.7 


0.5794043665 


9.1891232 


0.8 


0.6169244095 


15.8955019 


0.9 


0.6956250340 


47.7590095 


0.95 


0.7788797565 


166.7543222 


0.99 


0.9200729380 


3906.8991705 


0.991 


0.9259539055 


4823.955200 


0.992 


0.9320895460 


6107.378535 


0.993 


0.9385115510 


7981.5666867 


0.994 


0.9452603485 


10873.01788 


0.995 


0.9523894295 


15675.49526 


0.996 


0.9599731720 


24531.28407 


0.997 


0.9681226815 


43701.09808 


0.998 


0.9770238615 


98596.55052 


0.999 


0.9870615640 


395877.877549 


0.9999 


0.9982370395 0.3980357873 x 10^ 


1.7 


1.043246242 


0.6064576241 x 10^ 


3.0 


0.5120894280 


10410.51705 


5.0 


0.5023884065 


6235.986909 
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mc=0.51 70643255, L^=2.8743397865, Y=0.0 




FIG. 8: The potential V(R) for 7 = with some values of L near the critical point L = Lc- The 
curves from top to botton represent Li to L4, respectively. 

see Figs. 19-24. In particular, as 7 — > 1, we find that —>■ 00. This can be seen from Table 
III and Figs. 12-18. 

In the cases where we can have the two types of stable gravastars, it is also possible to 
find configurations where black holes are formed. This shows clearly that, even gravastars 
indeed exist, they do not exclude the existence of black holes. 



III. CONCLUSIONS 



In this paper, we have generalized our previous work on the problem of stable gravastars 
by constructing dynamical three-layer VW models jsj, which consists of an internal de 
Sitter space, a dynamical infinitely thin shell of a perfect fluid with the equation of state 
p = (1 — 7)0", and an external Schwarzschild spacetime. We have shown explicitly that 
the final output can be a black hole, a "bounded excursion" gravastar, a stable gravastar, 
a Minkowski, or a de Sitter spacetime, depending on the total mass m of the system, the 
cosmological constant A, and the initial position Rq of the dynamical thin shell. All these 
possibihties have non-zero measurements in the phase space of m, 0), 7(< 1) and i?o, 
although the region of gravastars is very small in comparison with that of black holes. When 
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0.4 
0.3 
0.2 
0.1 



£ 



-0. 

-0.2- 
-0.3- 
-Oj 




m^=0.500 

m3=0.530 
m4=0.550 



1.5 
R 



2.5 



FIG. 9: The potential V(R) for 7 = with some values of m near the critical point m = rric- The 
curves from top to botton represent mi to m^, respectively. 




FIG. 10: The potential V(R) for 7 = 0.4 with some values of L near the critical point L = Lc- The 
curves from top to botton represent Li to L4, respectively. 



16 



0.4 

0.3 

0.2 

0.1 

-0. 
-0.2- 
-0^ 



mc=0.5350989795, L^=4.3336020545, y=0.4 



mi=0.510 

m3=0.545 
m4=0.560 



1.5 



2 
R 



2.5 



3.5 



FIG. 11: The potential V(R) for 7 = 0.4 with some values of m near the critical point m = rric- 
The curves from top to botton represent mi to 771,4, respectively. 

7 > 1 even with A 7^ 0, only black holes are found. Therefore, although the existence of 
gravastars cannot be completely excluded in these dynamical models, our results show that, 
even if gravastars indeed exist, they do not exclude the existence of black holes. 
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Padmanabhan, Proc. of the 29th Int. Cosmic Ray Conf. 10, 47 (2005); E.J. Copeland, M. 
Sami, and S. Tsujikawa, ^^Dynamics of dark energy" |arXiv:hep-th/0603057 (2006); E.W. 



Kolb, Cosmology and the Unexpected," arXiv:0709.3102 : E. Linder, Mapping the cosmolog- 



20 



.x10 




1000 2000 3000 4000 5000 6000 7000 8000 

R 



FIG. 18: The potential V(R) for 7 = 0.9999 near its minimal point. 
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FIG. 19: The potential V(R) for 7 = 1.7. The curves from top to botton represent mi to m^, 
respectively. 
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FIG. 20: The potential V(R) for 7 = 1.7. The curves from bottom to top represent L4 to Li, 
respectively, where the curves Li and L2 coincide. Any potential curve where L > Lc will coincide 
with the potential curve where L = Lc 
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FIG. 21: The potential V(R) for 7 = 3.0. The curves from top to botton represent mi to m^, 
respectively. 
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FIG. 22: The potential V(R) for 7 = 3.0. The curves from bottom to top represent L4 to Li, 
respectively, where the curves Li and L2 coincide. Any potential curve where L > Lc will coincide 
with the potential curve where L = Lc- 
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FIG. 23: The potential V(R) for 7 = 5.0. The curves from top to botton represent mi to m^, 
respectively. 
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FIG. 24: The potential V(R) for 7 = 5.0. The curves from bottom to top represent L4 to Li, 
respectively, where the curves Li and L2 coincide. Any potential curve where L > Lc will coincide 
with the potential curve where L = Lc- 



